function shell_report_minloss(OO_minloss, ...
     D, K, N, D_enroll_st,  pr_y, alpha, B1_c, beta_c, T1, B2_c,...
     low_multiplier, high_multiplier, export_tex, folder_tex_cf, resource_y1, T2)

ssolve_market_all = OO_minloss(1).solve_market_all_cf; % N.ist x 1
ssolve_market     = OO_minloss(1).solve_market_cf; % N.st x 1

% check convergence under three scenarios
nnon_converge_all_rsr2000 = 0; % N.ist x 1
nnon_converge_rsr2000 = 0; % N.st x 1
for ss=1:3
    nnon_converge_all_rsr2000 = nnon_converge_all_rsr2000 + OO_minloss(ss).non_converge_all_cf;
    nnon_converge_rsr2000 = nnon_converge_rsr2000 + OO_minloss(ss).non_converge_cf;
end

ss_ind = 3;

p1_eq            = OO_minloss(ss_ind).p1_cf; % N.ist x 1
p2_eq            = OO_minloss(ss_ind).p2_cf; % N.ist x KM
n_eq             = OO_minloss(ss_ind).n_cf;  %N.st x 1
profit_eq        = OO_minloss(ss_ind).profit_cf; % N.ist x 1
non_converge_all = OO_minloss(ss_ind).non_converge_all_cf; % N.ist x 1
solve_market_all = OO_minloss(ss_ind).solve_market_all_cf; % N.ist x 1
fsolve_all       = OO_minloss(ss_ind).fsolve_all_cf; % N.ist x 1
EV_y             = OO_minloss(ss_ind).EV_y_cf; % N.st x N_y
enrollment_eq    = OO_minloss(ss_ind).enrollment_cf; % N.st x 1
enrollment_m_eq  = OO_minloss(ss_ind).enrollment_m_cf; % N.st x 1
enrollment_f_eq  = OO_minloss(ss_ind).enrollment_f_cf; % N.st x 1
non_converge     = OO_minloss(ss_ind).non_converge_cf; % N.st x 1
solve_market     = OO_minloss(ss_ind).solve_market_cf; % N.st x 1
rs_share_eq      = OO_minloss(ss_ind).rs_share_cf;  % N.ist x 1

[fit_off_st, fit_off_ist] = report_baseline(enrollment_eq, non_converge, solve_market,...
    D, K, N, D_enroll_st);

ss_ind = 3; 
ns = 2; % low regulatory cost & high regulatory cost will be considered
RS_rsr2000 = report_outcome(OO_minloss, ns, ss_ind, ...
    ssolve_market_all, nnon_converge_all_rsr2000, ssolve_market, nnon_converge_rsr2000,...
    D, K, N, D_enroll_st, pr_y, alpha, B1_c, beta_c, T1, B2_c,...
    fit_off_st, fit_off_ist, resource_y1, T2);

if export_tex==1
    file_tex  = fullfile(folder_tex_cf,'cf_minloss.tex');
    FID = fopen(file_tex, 'w');
    fprintf(FID, '\\begin{table} \n');
    fprintf(FID, '\\begin{center} \n');
    fprintf(FID, '\\scalebox{1.0}{ \n');
    fprintf(FID, '\\begin{tabular}{lcc} \\hline \n');
    fprintf(FID, ' & (1)  & (2) \\\\ \n');
    fprintf(FID, ' & Less strict & More strict \\\\ \n');
    fprintf(FID, '\\hline \n');
    fprintf(FID, 'Initial rate & %10.2f   & %10.2f \\\\ \n', RS_rsr2000.p1_vec);
    fprintf(FID, 'Revised rate & %10.2f   & %10.2f \\\\ \n', RS_rsr2000.p2_vec);
    fprintf(FID, '\\hline \n');
    fprintf(FID, '\\end{tabular} \n');
    fprintf(FID, '} \n');
    fprintf(FID, '\\end{center} \n');
    fprintf(FID, '\\caption{Price dynamics under counterfactual initial rate regulation \\newline \\small {\\emph{Notes}: The table reports percent changes in outcomes relative to the baseline economy as we vary the initial rate regulatory cost. In column (1), we decrease the initial rate regulatory cost to a half of the baseline estimate. In column (2), we double the cost.}} \n');
    fprintf(FID, '\\label{table_cf_minloss} \n');
    fprintf(FID, '\\end{table} \n');

end

